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Background: The underlying molecular mechanisms of PTSD are largely unknown. Dis- 
tinct expression signatures for PTSD have been found, in particular for immune activation 
transcripts. DNA methylation may be significant in the pathophysiology of PTSD, since the 
process is intrinsically linked to gene expression. We evaluated temporal changes in DNA 
methylation in select promoter regions of immune system-related genes in U.S. military 
service members with a PTSD diagnosis, pre- and post-diagnosis, and in controls. 

Methods: Cases (n = 75) had a post-deployment diagnosis of PTSD in their medical record. 
Controls (n = 75) were randomly selected service members with no PTSD diagnosis. DNA 
was extracted from pre- and post-deployment sera. DNA methylation (%5-mC) was quan- 
tified at specific CpG sites in promoter regions of insulin-like growth factor 2 (IGF2), long 
non-coding RNA transcript H19, interleukin-8 (IL8), IL16, and IL18 via pyrosequencing. We 
used multivariate analysis of variance and generalized linear models to calculate adjusted 
means (adjusted for age, gender, and race) to make temporal comparisons of %5-mC for 
cases (pre- to post-deployment) versus controls (pre- to post-deployment). 

Results: There were significant differences in the change of %5-mC pre- to post- 
deployment between cases and controls for H19 (cases: +0.57%, controls: -1.97%; 
p = 0.04) and IL18 (cases: +1.39%, controls: -3.83%; p = 0.01). For H19 the difference 
was driven by a significant reduction in %5-mC among controls; for I LI 8 the difference 
was driven by both a reduction in %5-mC among controls and an increase in %5-mC 
among cases. Stratified analyses revealed more pronounced differences in the adjusted 
means of pre-post H19 and I LI 8 methylation differences for cases versus controls among 
older service members, males, service members of white race, and those with shorter 
deployments (6-12 months). 

Conclusion: In the study of deployed personnel, those who did not develop PTSD had 
reduced %5-mC levels of H19 and I LI 8 after deployment, while those who did develop 
PTSD had increased levels of I LI 8. Additionally, pre-deployment the people who later 
became cases had lower levels of I LI 8 %5-mC compared with controls. These findings 
are preliminary and should be investigated in larger studies. 

Keywords: post-traumatic stress disorder, DNA methylation, cytokines, inflammation, promoter region. Operation 
Iraqi Freedom, Operation Enduring Freedom, epigenetics 



INTRODUCTION 

Estimates of post-deployment PTSD prevalence in U.S. military 
service members who served in Afghanistan and Iraq are 12 and 
18%, respectively (Hoge et al., 2004); a recent study reported aver- 
ages of 10-20% (Thomas et al., 2010). Despite these high rates 
among military populations, as well as high rates in civilian pop- 
ulations [i.e., U.S. prevalence rates of 7-8% (Kessler et al, 1994, 
1997); lifetime prevalence of 9-14% (Kessler et al, 1995; Breslau, 



2001)], the underlying molecular mechanisms of PTSD are largely 
unknown. 

Given the sustained stress of the sympathetic nervous system 
and resulting hyper-arousal state in PTSD, it has been hypoth- 
esized that immune system functioning is affected. A series of 
recent studies has shown that PTSD is associated with alterations 
in immune function (Hoge et al., 2009) in particular increased lev- 
els of some canonical pro-inflammatory cytokines (Spivak et al., 
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1997; Maes et al, 1999; Sutherland et al, 2003; von Kanel et al, 
2007) and decreased levels of anti-inflammatory cytokines (Kawa- 
mura et al, 2001; von Kanel et al, 2007; Smith et al, 2011). 
Profiling using cDNA microarrays of peripheral blood during 
the triggering and development of PTSD in trauma survivors at 
the emergency room and 4 months later has found differential 
gene expression signatures in promoter regions of genes which 
distinguished PTSD patients (Segman et al., 2005; Zieker et al., 
2007). Distinct expression signatures were found in particular for 
transcripts involved in immune activation (Segman et al., 2005). 
Up-regulated genes included insulin-like growth factor 2 (IGF2) 
and pro-inflammatory/neuroprotective chemokine, IL8 (CXCL8); 
down-regulated genes included the pro-inflammatory pleiotropic 
cytokines, interleukin 16 (IL16), and interleukin 18 (IL18) (Seg- 
man et al, 2005; Zieker et al., 2007). Similar findings of differential 
expression of genes involved in immune cell function were found 
among cases of PTSD from the World Trade Center attack ver- 
sus controls (Sarapas et al., 2011) and in a study of Bosnian war 
refugees with PTSD versus controls (Nowotny et al, 2010). 

An epigenetic mechanism, DNA methylation may play a sig- 
nificant role in the pathophysiology of PTSD, since the process is 
intrinsically linked to gene regulation. Recently, epigenetic mech- 
anisms, including regulation of chromatin structure and DNA 
methylation, have been found to be regulators of gene transcrip- 
tion in the CNS (Zovkic and Sweatt, 2013). Methylation changes 
may completely silence a gene or may decrease or increase gene 
expression. It is hypothesized that epigenetic molecular mecha- 
nisms, especially DNA methylation/demethylation, may influence 
long-term behavioral change through active regulation of gene 
transcription in the CNS. A recent review proposed that epigenetic 
molecular mechanisms underlie the formation and stabilization 
of context- and cue-triggered fear conditioning based in the hip- 
pocampus and amygdala, a conclusion reached in a wide variety 
of studies using laboratory animals (Zovkic and Sweatt, 2013). 
In humans, there is a small, but growing body of literature that 
supports a potential link between DNA methylation of immune 
function genes and PTSD. A recent study applying methylation 
microarrays to assay CpG sites in peripheral blood from PTSD 
cases and controls found differential methylation in genes related 
to immune system functions (Uddin et al, 2010). In particular, 
this study found that IL8, a gene that regulates innate and adaptive 
immune system processes, was unmethylated among PTSD cases 
(Uddin et al, 2010). Another recent study in an African Amer- 
ican population found that psychosocial stress may alter global 
and gene-specific DNA methylation patterns potentially associated 
with peripheral immune dysregulation (Smith et al, 2011). 

In the present study, we focused on a set of growth fac- 
tors/transcripts and cytokines/chemokines previously implicated 
in PTSD: paternally imprinted IGF2, maternally imprinted H19 (a 
long, non-coding RNA transcript), IL8, IL16, and IL18 (Segman 
et al, 2005; Zieker et al, 2007; Nowotny et al, 2010; Uddin et al, 
2010; Sarapas et al., 2011; Smith et al, 2011). We evaluated pro- 
moter region methylation levels as percentage of 5-methyl cytosine 
(%5-mC) in DNA from serum of PTSD cases and controls, who 
were U.S. military soldiers who deployed to Afghanistan [Opera- 
tion Enduring Freedom (OEF)] or Iraq [Operation Iraqi Freedom 
(OIF)] between 2004 and 2006. Our study had the unique ability to 



investigate temporal changes in methylation patterns after deploy- 
ment [a proxy for exposure to a potentially traumatic event (PTE)]. 
The pre- and post-deployment samples collected in our study 
enable such evaluations of temporal changes in methylation pat- 
ters. Our objectives were to investigate DNA methylation patterns 
associated with PTSD. 

Since human studies of brain tissue are highly invasive, identify- 
ing a low-invasive biomarker of epigenetic patterns of PTSD would 
be of great clinical value. Most of the studies to date which have 
measured expression signatures or methylation patterns in PTSD 
or other psychiatric disorders have been carried out using periph- 
eral blood (Segman et al., 2005; Zieker et al., 2007; Nowotny et al., 
2010; Uddin et al, 2010; Sarapas et al, 2011; Smith et al, 2011). 
To our knowledge, serum has not yet been evaluated for gene- 
specific methylation patterns potentially associated with PTSD. 
Serum and cerebrospinal fluid (CSF) have been found to have 
good correlation with respect to cytokine gene expression, and in 
this exploratory study we evaluated DNA methylation from serum 
as a biomarker, without drawing links to other types of tissues 
or directly extrapolating its significance. Serum DNA methylation 
patterns may provide a surrogate indicator of differential response 
to stress and vulnerability to PTSD. Unlike previous studies, we 
had unique access to biologic samples prior to onset of disease 
and prior to deployment (a proxy for exposure to a PTE). Serum 
samples from military service members have been stored at the 
Department of Defense Serum Repository (DoDSR) since the 
early 1990's. We were able to access those stored sera and iden- 
tified samples drawn prior to and post-deployment for all cases 
and controls. 

MATERIALS AND METHODS 
STUDY POPULATION 

The study population has been described previously (Rusiecki 
et al., 2012). The target population was male and female U.S. 
Army and Marines service members serving their first OEF/OIF 
deployment between January 01, 2004 and December 31, 2006. 
Deployment length was between 6 and 18 months. Via query of 
medical records using the International Classification of Diseases, 
Ninth Revision (ICD-9) codes 290-320, we determined an absence 
of any mental health diagnoses dating back to at least 2 years 
prior to first OEF/OIF deployment for all cases and controls. In 
an effort to exclude possible confounding by other psychiatric 
illnesses, for which differential gene expression or methylation 
has been reported (Mill and Petronis, 2007; Shimabukuro et al., 
2007; Tamura et al., 2007; Kuratomi et al., 2008; Feng and Fan, 
2009; Iwamoto and Kato, 2009; Gavin and Sharma, 2010) post- 
deployment exclusion criteria for both cases and controls was ever 
having a health encounter for schizophrenia (ICD-9 code 295), 
bi-polar disorders (ICD-9 code 296), and manic phase bi-polar 
disorder (also ICD-9 code 296). 

The PTSD cases (« = 75) had existing pre- and post- 
deployment serum samples housed at the DoDSR, met all the 
criteria above, and had at least two outpatient records with a pri- 
mary diagnosis of chronic PTSD, based on ICD-9 Code 309.81 
in the first diagnostic position. The first outpatient diagnosis was 
between 4 and 12 months after return from first deployment. The 
second outpatient diagnosis was any time after that, but within 
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2 years of return from first deployment. Additional criteria for 
inclusion as a case to this study was having one serum sample 
(the pre-deployment sample) drawn within 12 months prior to 
first OEF/OIF deployment and one sample (the post-deployment 
sample) drawn within 6 months after return from first OEF/OIF 
deployment. Cases were randomly selected within each gender, 
such that 25 females and 50 males were included. We identified an 
appropriate active duty service member control group (n = 75), 
who were frequency matched from a stratified, random sample 
based on age (20-26, 27-35), gender, and race, but for whom there 
was never a diagnosis of PTSD (ICD-9 Code = 309.81) or Trau- 
matic Brain Injury (TBI) (ICD-9 Codes = 800.0-801.9, 803.0- 
804.9, or 850.0-854.1). Though cases in this study are presumed 
to have experienced a PTE by virtue of their PTSD diagnosis, the 
criteria of which requires exposure to a PTE, we have no data on 
formal assessment of trauma/PTE for either cases or controls. As 
such, this study uses deployment as a proxy for exposure to a PTE 
for both cases and controls. Because of Department of Defense 
policy regarding the DoDSR samples, we were unable to make con- 
tact with the study subjects for whom we had serum and medical 
encounter data. 

SAMPLE PREPARATION AND LABORATORY METHODS 
DNA extraction 

For each PTSD case and control, The Armed Forces Health Surveil- 
lance Center identified a pre-deployment and a post-deployment 
serum sample and authorized release of up to 0.5 ml of serum per 
sample. These samples had been maintained at —30 °C contin- 
uously at the DoDSR. DNA was extracted using ChargeSwitch 
gDNA 0.2-1 ml serum kit from (Invitrogen Carlsbad, CA, USA) 
and quantified via Qubit 8 dsDNA HS Assay Kit using a Qubit 
fiuorometer (Invitrogen, Carlsbad, CA, USA). 

Quantification of DNA methylation 

DNA methylation was quantified via bisulfite treatment, PCR, and 
pyrosequencing. DNA was bisulfite treated using the Zymo DNA 
Methylation Kit (Zymo research, Orange, CA, USA). All assays 
used in this study were validated by PCR bias testing and sensitiv- 
ity testing. This was done to assess the lowest input gDNA needed 
for bisulfite modification and the lowest input bisulfite converted 
DNA for PCR for each assay. For H19, IL8, IL16, and IL18, a mini- 
mum of 40 ng of DNA was used for bisulfite conversion. For IGF2, 
a minimum of 10 ng DNA was used for bisulfite conversion. The 
total DNA available for each sample was used as input gDNA for 
bisulfite conversion. Carrier RNA from Qiagen was used during 
the bisulfite modification and purification to minimize the loss of 
DNA during the procedure. Bisulfite treated DNA was eluted in 
20 |xl volume and 1 |xl of it was used for each PCR. The PCR (45 
cycles) was performed with one of the PCR primers biotinylated to 
convert the PCR product to single-stranded DNA templates. The 
PCR products (each 10 |xl) were sequenced by Pyrosequencing 
PSQ96 HS System (Qiagen Pyrosequencing) following the manu- 
facturer's instructions (Qiagen Pyrosequencing). Pyrosequencing 
is a real-time sequencing-based on mutation analysis or methy- 
lation analysis technology. The methylation status of each CpG 
locus was analyzed individually as a T/C SNP using QCpG software 
(Qiagen Pyrosequencing). 



The loci of specific CpGs measured in each promoter region 
are shown in Figure 1 and the rationale for selecting these spe- 
cific loci was that they are generally located in the shore/vicinity 
of transcriptional start site. Loci measured were in the promoter 
regions of H19 [-1964, -1946, -1927, -1919 from the transcrip- 
tional H19 start site (TSS)], IGF2 [TSS-479, -476, -460, -361, -341, 
-322; note, the investigated regions for IGF2 covered the promoter 
of IGF2, p2, and the binding region for enhancer-blocking ele- 
ment CCCTC-binding factor (CTCF) upstream of the H19 start 
site], IL8 (TSS-116, -106, -31), IL16 (TSS-159, -139, -93, -79), and 
IL18 (TSS-158, -108, -86, -49, -33). The percentage of methylation 
(%5-mC) was expressed as 5-mC divided by the sum of methylated 
and unmethylated cytosine. This is interpreted as the percentage 
of cytosines at a given CpG site (or position) which is methy- 
lated. Percent 5-mC was measured at each CpG site, hereafter 
referred to as position, and a mean %5-mC was calculated across 
all positions measured in a promoter region. Starting from bisulfite 
modification, four controls [low, medium, high methylated DNA 
(EpigenDx, Inc.), and a no DNA template] were included in every 
pyrosequencing run to ensure specificity of PCR amplification, 
and success of pyrosequencing reactions. 

DNA from cases and controls were randomly arranged on the 
plates for methylation quantification analysis. For 10% of the sam- 
ples we included duplicates to which laboratory personnel were 
blinded, for the purpose of quality control (QC). Coefficients of 
variation as well as intraclass correlation coefficients (ICCs) were 
calculated for the means of each position in H19 (CV = 0.05, 
ICC = 0.75; £<0.01), IL8 (CV = 0.01, ICC = 0.72; p<0.01), 
IL16 (CV = 0.05, ICC = 0.68; p<0.01), and IL18 (CV = 0.02, 
ICC = 0.89; p < 0.01). The duplicate samples were run on a sepa- 
rate plate from their counterparts, and some of the variation could 
be due to a plate effect. 

Approximately 30% of the samples (N = 86) had DNA yields 
(e.g., <40ng) insufficient for bisulfite conversion of HI 9, IL8, 
IL16, and IL18. Approximately 5% of the samples (N = 15) had 
DNA yields (e.g., <10ng) insufficient for bisulfite conversion of 
IGF2. These samples were excluded from our study. For a small 
percentage of the samples for which bisulfite conversion was com- 
plete, pyrosequencing signals were too low to be reliable (0% for 
H19, 9% for IGF2, 0.05% for IL8, 8% for IL16, and 2% for IL18), 
thus we did not include them in our final analyses. It is unclear 
why signals were low for those particular samples, however, low 
signals were not associated with DNA yield. Distribution of sam- 
ples either not measured or measured but with low signals did not 
differ by pre- and post-deployment or case/control status, with the 
exception of IGF2 pre-deployment case/control status. However, 
this is likely because there were only 16 with low signals for that 
assay, thus lending to the instability of those proportions. 

Statistical methods 

Our statistical analyses considered %5-mC in the promoter and/or 
known methylation region of each gene at each position (CpG 
site) measured as well as the mean of all the positions mea- 
sured. We compared the change in cases' %5-mC from pre- to 
post-deployment, to the change in controls' %5-mC from pre- 
to post-deployment using multivariate ANOVA. We present the 
adjusted mean of pre-post methylation difference for cases versus 
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FIGURE 1 | Loci of specific CpGs measured in IGF2, H19, IL8, IL16, and IL18 



controls, adjusted for age, gender, and race. To further investi- 
gate the adjusted means and differences between cases and con- 
trols, pre-deployment and post-deployment, we ran generalized 
linear models (GLMs) Using GLMs, we carried out case-case 
and control-control (i.e., pre- to post-deployment) comparisons, 
as well as case-control comparisons, pre-deployment, and post- 
deployment. We carried out the same analyses, stratifying by age, 
gender, race, and length of deployment, separately, to investigate 
any potential patterns within specific groups. SAS Version 9.2 was 
used to carry out all statistical analyses. The SAS procedures GLM 
and ANOVA were used to run our final models (SAS Institute Inc., 
2008). 

This study was approved by the Institutional Review Board at 
the Uniformed Services University of the Health Sciences. 

RESULTS 

Baseline characteristics of the study population are in Table 1. 
Although the study was designed with 75 cases and 75 controls, 
we excluded 2 study subjects (1 case/1 control) because their 
DNA yield was insufficient for DNA methylation quantification 



for any of the loci measured in both pre- and post-deployment 
time frames. The remaining 148 study subjects did not neces- 
sarily have complete data for all loci measured in both pre- and 
post-deployment time frames, but they did have DNA methyla- 
tion quantified for at least one locus in at least one time frame. 
The study population did not differ by case-control status for age, 
gender, and race because of the selection and frequency match- 
ing criteria. Approximately 73% of cases and 76% of controls had 
deployments of 6 months to less than 12 months, while 27% of 
cases and 24% of controls had deployments of 12-18 months. 
These small differences were not statistically significant. The num- 
ber of days between end of deployment and post-deployment 
serum draw ranged from 1 to 170 (mean = 22.45; SD = 38.84; 
median = 7), and the number of days between pre-deployment 
serum draw and start of deployment ranged from 58 to 358 
(mean = 87.46; SD = 84.49; median = 58.5). There were no dif- 
ferences between cases and controls for these time intervals, nor 
were there differences by occupational code or rank/rate. There 
were four cases who also had a post-deployment diagnosis of 
TBI (defined as an ICD-9 code of 800.0-801.9, 803.0-804.9, or 
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Table 1 | Baseline characteristics of population. 
Characteristic 



Cases (« = 74) 



N 



Controls (« = 74) 



N 



p-Value 





Younger (20-23 years) 


40 


54.1 


38 


51.4 


0.74 


Older (24-35 years) 


34 


46.0 


36 


48.7 






Male 


49 


66.2 


49 


66.2 


1.00 


Female 


25 


33.8 


25 


33.8 




White 


59 


79.7 


60 


81.1 


0.84 


Black 


15 


20.3 


14 


18.9 




Short (6 to <12 months) 


54 


73.0 


56 


75.7 


0.71 


Long (12-18 months) 


20 


27.0 


18 


24.3 




< 7 days 


41 


55.4 


38 


51.4 


0.62 


> 7 days 


33 


44.6 


48 


48.6 




< 90 days 


52 


70.3 


45 


60.8 


0.23 


> 90 days 


22 


29.7 


29 


39.2 




OCCUPATIONAL CODING* 


Less likely involved in combat 


19 


38.8 


20 


38.5 


0.97 


More likely involved in combat 


30 


61.2 


32 


61.5 




Missing 


25 




22 








Junior enlisted 


57 


77.0 


52 


70.3 


0.35 


Middle-grade/senior enlisted and officer 


17 


23.0 


22 


29.7 





'Based on military occupational specialty POC coding as either less likely to be involved in combat or more likely to be involved in combat. 
"Junior enlisted included E-2 through E-4; middle/senior enlisted and officer included E-5 through E-8, O- 1 through 0-3, and W-2. 



850.0-854.1) (data not shown). We carried out all analyses includ- 
ing and excluding those four cases to evaluate potential confound- 
ing by TBI and found that results were very similar. We, therefore, 
present all results for analyses including those four PTSD cases 
with concurrent TBI. 

Table 2 presents the results of a multivariate ANOVA com- 
paring adjusted means of pre-post methylation differences for 
cases versus controls. These results show that, accounting for the 
exposure of deployment (i.e., proxy for PTE), there were signif- 
icant reductions in %5-mC for controls, compared with cases at 
H19 (mean of positions 1-4: cases: +0.57%, controls: —1.97%; 
£ = 0.04) and IL18 (mean of positions 1-5: cases: +1.39%, con- 
trols -3.83%; p = 0.01). Position 3 in H19 and positions 4 and 5 
in IL18 had statistically significant differences between cases and 
controls, and though differences in the change of methylation were 
not significant at other positions for H19 and IL18, they all showed 
similar patterns. 

Table SI in Supplementary Material provides additional detail 
on the change in methylation across deployment for cases and con- 
trols; paired differences in adjusted mean levels for cases pre- and 
post-deployment and controls pre- and post-deployment, deter- 
mined via GLM, are presented. The differences between cases and 



controls across deployment for HI 9 were driven primarily by the 
significant reduction in %5-mC among controls from pre- to post- 
deployment (the levels for cases were unchanged across deploy- 
ment for H19). The pattern of decreased H19 methylation among 
controls from pre- to post-deployment is evident; statistically sig- 
nificant decreases were found at positions 1, 3, and 4, and for 
the mean of positions 1-4. The differences for IL18 were dri- 
ven by both the decrease in methylation across deployment for 
controls as well as the increase in methylation across deployment 
for cases. There was a pattern of decreased methylation among 
controls from pre- to post-deployment, though a statistically sig- 
nificant decrease was found for position 5 only. There was also 
a pattern of increased IL18 methylation for cases pre- to post- 
deployment, though none of the comparisons was statistically 
significant. 

Table S2 in Supplementary Material, which presents case- 
control comparisons for both pre- and post-deployment, further 
illustrates that for IL18, %5-mC was lower in cases compared with 
controls, pre-deployment (p < 0.05 for positions 2, 4, and 5, and 
the mean of the five positions), while the pattern reversed post- 
deployment, such that levels were non-statistically significantly 
higher in cases compared with controls. Note that the term "case" 
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Table 2 | Multivariate ANOVA comparing adjusted* mean of pre-post 
methylation difference between cases and controls. 



Gene promoter 
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37 


0.57 


1.10 


0.04 




Control 
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0.99 
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0.09 




Control 


45 


-1.92 


1.16 
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Case 


37 


1.38 


1.41 


0.14 




Control 


45 


-1.01 


1.26 




Position 3;TSS-1926 


Case 


37 


0.65 


1.15 


0.01 




Control 


45 


-2.93 


1.03 
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Case 


37 
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0.31 




Control 
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1.29 
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Case 


52 
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0.40 




Control 
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Case 


52 


-12.14 


4.59 


0.67 




Control 


62 


-9.98 


4.32 




Position 5;TSS-341 


Case 


52 


-2.79 


3.97 


0.84 




Control 


62 


-3.71 


3.74 




Position 6;TSS-322 


Case 


52 


-2.90 


3.37 


0.95 




Control 


62 


-3.13 


3.17 




IL8 (mean of positions 1-3) 


Case 


36 


-0.17 


0.32 


0.25 




Control 


44 


-0.58 


0.28 




Position 1;TSS-116 


Case 


36 


-0.54 


0.54 


0.38 




Control 


44 


-1.06 


0.47 




Position 2;TSS-106 


Case 


36 


-0.12 


0.46 


0.36 




Control 


44 


-0.58 


0.40 




Position 3;TSS-31 


Case 


36 


0.15 


0.17 


0.20 




Control 


44 


-0.10 


0.15 




IL16 (mean of positions 1-4) 


Case 


31 


-4.23 


3.66 


0.88 




Control 


38 


-3.61 


3.66 




Position 1;TSS-159 


Case 


31 


-5.61 


4.80 


0.97 




Control 


38 


-5.44 


4.78 




Position 2;TSS-139 


Case 


31 


-2.94 


3.34 


0.89 




Control 


38 


-3.44 


3.33 




Position 3;TSS-93 


Case 


31 


-4.29 


6.58 


0.95 




Control 


38 


—4.76 


6.64 




Position 4; TSS-79 


Case 


31 


-5.00 


6.22 


0.85 




Control 


38 


-3.78 


6.27 




I LI S (mean of positions 1-5) 


Case 


37 


1.39 


1.85 


0.01 




Control 


45 


-3.83 


1.66 




Position 1;TSS-158 


Case 


37 


-2.60 


2.66 


0.50 




Control 


45 


-4.63 


2.40 




Position 2;TSS-108 


Case 


37 


2.31 


2.31 


0.03 




Control 


45 


-3.50 


2.08 





(Continued) 



A/pairs 1 Mean 2 SE p- Value 



Position 3;TSS-86 


Case 


37 


1.88 


2.47 


0.21 




Control 


45 


-1.66 


2.22 




Position 4; TSS-49 


Case 


37 


2.35 


1.93 


0.01 




Control 


45 


-3.43 


1.75 




Position 5;TSS-33 


Case 


37 


2.13 


2.23 


0.02 




Control 


45 


-4.03 


2.02 





'Adjusted for age, gender, race. 

' Pairs of samples, i. e., pre and post paired samples among the cases and among 
the controls - number of pairs differ for the different promoter regions measured 
because of sensitivity of each assay -i.e., some samples were dropped because 
of low pyrosequencing signals for a specific assay. 

2 Adjusted mean of the difference from pre- to post-deployment for cases and for 
controls. 

in the pre-deployment context refers to people who later became 
cases, post-deployment; during the pre-deployment time period, 
they were not yet cases. 

Also shown in Table S2 in Supplementary Material are a few 
other loci where significant differences in %5-mC were found 
between cases and controls pre-deployment, i.e., IGF2 position 
2 (%5-mC C a S es > %5-mC controls* p = 0.01) and IL8 position 2 
(%5-mC cases < %5-mC controls) P — 0.01). There were no signifi- 
cant differences in adjusted mean levels for the post-deployment 
case-control comparisons. 

Stratified analyses (Table 3) revealed more pronounced differ- 
ences in the adjusted means of pre-post H19 and IL18 methylation 
differences for cases versus controls among older (24-35 years) 
service members (H19 position 3 case/control: +0.13/— 3.66; 
p = 0.05), males (H19 position 3 case/control: +0.80/-2.99; 
p = 0.01; IL18 mean case/control: +0.47/— 5.38), service mem- 
bers of white race (HI 9 mean case/control: +0.72/— 2.70; p = 0.01; 
IL18 mean case/control: +2.82/-2.16; p = 0.04), and those 
with deployments between 6 and 12 months (HI 9 position 3 
case/control: +0.16/— 3.20; p = 0.02; IL18 mean case/control: 
+2.52/-4.65;p = 0.01). 

DISCUSSION 

Our study found that there were significant differences in the 
change of methylation (%5-mC) across deployment between cases 
and controls for H19 and IL18. While there was no change in %5- 
mC in H19 for cases from pre- to post-deployment, there was a 
significant decrease among controls driving the difference in %5- 
mC change between the cases and controls. For IL18 there was 
a significant difference between the %5-mC decrease in controls 
and the %5-mC increase in cases, from pre- to post-deployment. 
These findings were more pronounced in both HI 9 and IL18 in 
younger service members, males, service members of white race, 
and those with shorter deployment length. There were a few addi- 
tional statistically significant findings at specific CpG sites, such 
as pre-deployment case-control differences in IGF2 position 2 
and in IL8 position 2. Given the low signals from the IGF2 assay 



Wpairs Mean 2 SE p-Value Gene promoter 
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Table 3 | Multivariate ANOVA comparing adjusted* mean of pre-post methylation difference for cases versus controls, stratified by age, gender, 
race, and deployment length. 

Gene promoter/CpG site Subject Age: younger (20-23) and older (24-35) 



Younger Older 







N 


H 


p- Value 


N 


M- 


p-Value 


H19 position 3 


Case 


20 


0.71 


0.11 


17 


0.13 


0.05 




Cont 


25 


-2.47 




20 


-3.66 




H19 mean of positions 


Case 


20 


0.34 


0.36 


17 


0.58 


0.07 




Cont 


25 


-1.34 




20 


-3.00 




IL18 mean of positions 


Case 


20 


0.60 


0.22 


17 


1.38 


0.05 




Cont 


25 


-3.94 




20 


-3.07 












Gender: female and male 












Female 






Male 








N 


M- 


p- Value 


N 


V- 


p-Value 


H19 position 3 


Case 


9 


0.59 


0.24 


28 


0.80 


0.01 




Cont 


16 


-2.74 




29 


-2.99 




H19 mean of positions 


Case 


9 


0.98 


0.12 


28 


0.91 


0.14 




Cont 


16 


-2.84 




29 


-1.27 




IL18 mean of positions 


Case 


9 


1.12 


0.35 


28 


0.47 


0.03 




Cont 


16 


-1.80 




29 


-5.38 












Race: black and white 












Black 






White 








N 


H 


p- Value 


N 


V- 


p-Value 


H19 position 3 


Case 


6 


-0.65 


0.43 


31 


1.64 


0.01 




Cont 


8 


-3.07 




37 


-2.30 




H19 mean of positions 


Case 


6 


-0.95 


0.73 


31 


0.72 


0.01 




Cont 


8 


0.14 




37 


-2.70 




IL18 mean of positions 


Case 


6 


-0.31 


0.21 


31 


2.82 


0.04 




Cont 


8 


-5.53 




37 


-2.16 










Deployment length: long (6 to 


<12) and short (12-18) 










Long Depl. 






Short Depl. 








N 


H 


p- Value 


N 


M- 


p-Value 


H19 position 3 


Case 


8 


4.20 


0.42 


29 


0.16 


0.02 




Cont 


11 


0.75 




34 


-3.20 




H19 mean of positions 


Case 


8 


3.56 


0.38 


29 


0.12 


0.12 




Cont 


11 


0.68 




34 


-2.14 




IL18 mean of positions 


Case 


8 


-2.08 


0.93 


29 


2.52 


0.01 




Cont 


11 


-2.46 




34 


-4.65 





"For analyses stratified on age, adjusted for gender and race; for analyses stratified on gender, adjusted for age and race; for analyses stratified on race, adjusted for 
age and gender; for analyses stratified on deployment length, adjusted for age, gender, and race. 
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being differentially distributed between pre-deployment cases and 
controls, we are cautious to interpret the results for IGF2. Whether 
these patterns represent markers of vulnerability or resiliency is 
speculative. Psychological stress incurred during deployment may 
be associated with a methylation response in these genes, but given 
the lack of PTE data for controls, this study is not able to address 
whether the response is protective of PTSD or puts one at risk 
for PTSD. 

The most consistent result from our study, which is detailed in 
Table S2 in Supplementary Material, is that cases pre-deployment 
had significantly lower %5-mC in the IL18 promoter region than 
controls pre-deployment. This difference was found at positions 
2, 4, and 5, and the mean of all five positions; p-values approached 
significance for positions 1 (0.08) and 3 (0.06). There was a rever- 
sal of this pattern post-deployment, in that cases had higher levels 
at all positions, compared to controls, though differences were 
not statistically significant. Given that DNA methylation is often 
inversely correlated with gene expression, in particular for immune 
system-related genes (Oliveira et al, 2009) these results for IL18 are 
generally consistent with findings from a previous study based on 
cDNA microarray investigation, which reported down-regulation 
of IL18 in PTSD cases (Segman et al, 2005). Interestingly, pro- 
inflammatory cytokine IL18, a key mediator of inflammation and 
immune response inducing interferon-y (IFN-y), is expressed in 
the brain and plays a significant role in a number of neuropatho- 
logical disorders (Bossu et al, 2010; Anderson and Rodriguez, 
2011). IL18 has been recently shown to attenuate breaks in the 
blood-brain barrier via IFN-y independent pathway, suggesting 
its potential neuroprotective role (Jung et al, 2012). Our results of 
increased IL18 methylation in cases would imply down-regulation 
of the gene, which is consistent with previous findings (Zieker et al., 
2007). 

Our findings also suggest a potential role of H19 in stress 
response (23). The imprinted IGF2-H19 locus (CTCF binding 
cite) encodes the growth promoting hormone IGF2 and a long 
non-coding RNA H19. The maternally inherited locus transcribes 
untranslated H19 RNA, which plays the role of a trans-regulator 
of the imprinted gene network controlling embryonic growth and 
was recently shown to be processed into micro-RNA (miRNA)-675 
(Cunningham et al., 2010; Keniry et al., 2012). Expression of both 
genes is partially regulated by a 3' distal enhancer, and methylation 
of the CTCF binding site on the paternal chromosome prevents 
binding of the CTCF insulator allowing for activation of the IGF2 
promoter. Lack of methylation of the CTCF binding site on the 
maternal chromosome, however, prevents IGF2 promoter activa- 
tion (Fu et al, 2004; Cunningham et al., 2010; Keniry et al, 2012). 
Although a previous study reported that IGF2 was up-regulated 
in the whole blood (Zieker et al., 2007) IGF2-mediated signal- 
ing has been implicated in fear extinction and was proposed as a 
therapeutic venue to attenuate excessive fear memory potentially 
via promotion of neuronal survival and/or self-renewal of neural 
stem cells through its interaction with the insulin receptor (Agis- 
Balboa et al, 201 1; Ziegler et al, 2012). This could suggest that it is 
not only IGF2 that may play a direct role in PTSD fear extinction, 
but H19 may be involved in the regulation of IGF2. Also, H19 or 
its miRNA products and inhibition of their putative targets may 
contribute to the stress response (23). 



Previously, we investigated DNA methylation in repetitive ele- 
ments, LINE-1 and Alu, in this population and found differ- 
ential methylation between pre- and post-deployment controls 
and in cases versus controls, post-deployment. In the current 
study, we selected specific cytokine promoter regions of interest 
based on previous findings in the literature of differential gene 
expression or methylation of those cytokines between PTSD cases 
and controls (Segman et al, 2005; Zieker et al., 2007; Nowotny 
et al, 2010; Uddin et al, 2010; Sarapas et al, 2011; Smith et al, 
2011). Increasing evidence for contribution by the chemokine 
network (Semple et al., 2010) and the known role of IGF2 (Torres- 
Aleman, 2010) in brain and vasculature development indicates 
it as an important target due to possible pre-morbid conditions. 
The neuro-chemokine IL8 has been implicated in brain devel- 
opment, neuroinfiammation, and synaptic transmission (Rostene 
et al., 2011). Cytokine IL16, a chemo-attractant for immune 
cells expressing surface CD4 molecules (e.g., T-cells) (Cruikshank 
and Little, 2008) targets a number of neurophysiological mem- 
brane proteins, such as immediate-early gene c-Fos, potassium 
and calcium channels, N-Methyl-D-aspartate (NMDA) -receptor 
subunits, and neuronal phosphatases through its PDZ domain 
and exerts its autocrine function in the brain with a distinct 
role in neuropoiesis (Kurschner and Yuzaki, 1999; Schwab et al, 
2001; Bannert et al., 2003). Chemokine IL8 and pro-inflammatory 
cytokine IL18 are known to induce inflammatory mediator inter- 
feron y and have been found to play neuroprotective roles in 
models of neurodegenerative conditions (Ryu et al, 2010; Semple 
etal.,2010). 

Only a handful of recent studies have evaluated gene-specific 
DNA methylation and PTSD (Uddin et al., 2010, 2011; Ressler 
et al., 201 1; Trollope et al., 2012) to our knowledge, and these stud- 
ies have used only post-PTE blood samples. A cross sectional study 
of PTSD-affected and -unaffected individuals enrolled in a longi- 
tudinal study which investigated methylation and immune func- 
tion profiles in DNA derived from whole blood in a methylation 
microarray reported that immune system functions were signif- 
icantly overrepresented among the genes uniquely unmethylated 
in those with PTSD (Uddin et al, 2010). This signature included 
IL8, a gene which regulates innate and adaptive immune system 
processes (Uddin et al., 2010). There is emerging evidence that 
chemokines like IL8 and its receptors not only regulate immune 
cell infiltration but also contribute to a variety of physiological 
functions underlying neurotransmission, neuro-protection, and 
neurogenesis (Rostene et al, 2011) and therefore may exert a 
neuro-modulatory role by changing the balance of neurotrans- 
mitters and affecting fear memory (Herry et al, 2008; Joels et al., 
2008). However, a peripheral blood mononuclear cell (PBMC) 
gene expression profile carried out in trauma survivors found that 
IL8 was under-expressed in PTSD cases, which would indicate 
the potential for hypermethylation of IL8 in these cases (Segman 
et al, 2005). Another recent human study found that persons who 
experienced traumatic events were at increased risk for PTSD, but 
only those with lower methylation levels of a serotonin transporter 
gene, SLC6A4. At higher methylation levels, individuals with more 
traumatic events were protected from PTSD (Koenen et al., 201 1). 
Our study did not find a post-deployment case-control difference 
in this gene. We did find a pre-deployment case-control difference 
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in one position in IL8 (position 2) in that people who later became 
cases (pre-deployment cases) had lower levels of IL8 than did 
controls. This could indicate susceptibility to a PTE or to the devel- 
opment of PTSD, however we found it in only one position in the 
promoter region. 

Serum has not previously been evaluated as a biomarker for 
DNA methylation patterns in cytokine gene promoter regions 
associated with PTSD. Most of the human studies which have 
investigated DNA methylation in PTSD have utilized predom- 
inantly whole blood-derived DNA (Uddin et al, 2010, 2011; 
Koenen et al, 2011; Ressler et al, 2011; Trollope et al, 2012). 
How correlated DNA methylation levels are in serum and whole 
blood to brain and other CNS tissues is not clear. However, it 
has been reported that some cytokines, such as IL18 can cross 
the blood-brain barrier (Alboni et al., 2010; Jung et al, 2012). 
The origin of circulating nucleic acids such as cell-free circulating 
DNA (cfcDNA) has been hypothesized to stem from necrosis or 
apoptosis (Gormally et al, 2007) though others report the pos- 
sibility of an active release from cells (Anker et al, 1999; Stroun 
et al, 2001; Goebel et al, 2005; Rhodes et al, 2006). Excessive 
stress may induce DNA damage in the form of oxidized nucleo- 
sides, strand breaks, apoptosis, and necrosis and may be a source 
of cfcDNA in serum (Atamaniuk et al., 201 1). Additionally, in the 
CNS, severe life stress leads to oxidative stress, and the accumu- 
lation of reactive oxygen species (ROS) is known to increase the 
susceptibility of brain tissue to damage (Schiavone et al, 2013). 
ROS, in turn, may also trigger molecular cascades leading to the 
blood-brain barrier permeability and neuronal death (Gu et al, 
2011). Compared with cultured cells, clinical specimens, such as 
whole blood, serum, and even brain tissue and other CNS tissues, 
contain a heterogeneous mixture of cell types, each contributing its 
own unique methylation profile to the final analysis. Indeed DNA 
methylation patterns have been found to differ globally or locally, 
depending on the brain region/sub-region of focus (Iwamoto and 
Kato, 2009). One study demonstrated that the human cerebral 
cortex has an epigenetic signature distinct from the cerebellum 
(Ladd-Acosta et al, 2007) while others have shown that global 
DNA methylation differs across sub-regions of the rat hippocam- 
pus (Miller and Sweatt, 2007; Brown et al., 2008). These results 
taken collectively indicate that even different cortical cells are 
likely to have distinct epigenomic patterns. Serum DNA methy- 
lation patterns may provide a surrogate indicator of differential 
response to stress and PTSD, and our intent was to evaluate it as 
a non-invasive biomarker without drawing links to other types 
of tissues or directly extrapolating its significance. We evaluated 
difference in storage time of serum as a potential factor which 
may have influenced DNA methylation patterns. All samples for 
this study were obtained from the DoD Serum Repository in June 
2009, and DNA was extracted from them soon thereafter. The 
time frame of serum being collected from our study subjects, 
based on the study sample criteria, was approximately January 
2003 through June 2007. We found that cases and controls did not 
differ, either pre-deployment or post-deployment with respect to 
storage time. The mean storage time for cases pre-deployment was 
4.5 years, for controls it was 4.8 years (p = 0.15). The mean storage 
time for cases post-deployment was 3.3 years, for controls it was 
3.7 years (p = 0.14). As we expected, the storage time did differ 



by pre- and post-deployment time period. Thus for cases pre- 
deployment the mean was 4.5 years, while post-deployment it was 
3.4 years (p = < 0.01). Likewise, the mean storage time for con- 
trols pre-deployment was 4.9 years, while post-deployment it was 
3.7 (p < 0.01). Based on these differences, we adjusted all ourmod- 
els in which we carried out a pre/post comparison (i.e., Table 2; 
Table SI in Supplementary Material) by storage time. However, 
after adjusting our models for difference in storage time between 
pre- and post-deployment samples (for both cases and controls), 
there were no changes in our adjusted means and ^-values, so we 
opted for a more parsimonious model which did not include this 
variable. 

Limitations of this study include lack of detailed information 
on deployment exposures for both the cases and controls. Deploy- 
ment was used as a proxy for a PTE, and the exact timing of a 
PTE, if it occurred, is not known. There is also the issue that cases, 
by virtue of their PTSD diagnosis, were presumably exposed to a 
PTE, while we do not know the same for controls. In an attempt 
to investigate a possible bias from this disparity, we re-ran all our 
analyses utilizing a control group comprised only of those likely 
involved in combat during deployment (the group specified in 
Table 1 as "more likely involved in combat"). The methylation 
patterns we found did not differ from those using the complete 
control group (data not shown) . We also tried to control for poten- 
tial differences in deployment experience by ensuring that all cases 
and controls had not been previously deployed, that they were all 
active duty Army or Marines, and that they were deployed from 
6 to 18 months, but there is still potential for significant variation 
among all our subjects with respect to intensity of combat during 
deployment. We do not have any data on previous exposure to a 
PTE or underlying personality traits, which have been hypothe- 
sized to influence vulnerability to PTSD (van Zuiden et al., 201 1). 
Military personnel, prior to deployment, undergo a medical exam 
and are required to complete a pre-deployment health assessment 
survey (Department of Defense form 2795) to ensure medical 
readiness for deployment. This would have applied to all subjects 
in our study. 

The timing of sample collection was not standardized, so there 
was heterogeneity in the length of time between deployment and 
serum draw. Although we tried to minimize this interval, the tim- 
ing of each serum sample being added to the DoDSR depended 
on the timing of HIV testing for each military service member. 
We do not have data on other relevant exposures which may influ- 
ence DNA methylation, such as dietary factors (folate, vitamin Bn 
intake) (Fenech, 2001; Piyathilake and Johanning, 2002) smoking 
(Toh et al, 2010), and alcohol consumption (Seitz and Stickel, 
2007; Toh et al, 2010; Zhu et al., 2012). 

Ascertainment of PTSD via query of ICD-9 coded medical 
encounter is not ideal. Although we attempted to restrict the 
definition of PTSD to a scenario which would minimize misclassi- 
fication of disease, by requiring the ICD-9 code 309.81 be present 
in the first diagnostic position for two outpatient records spaced 
at a reasonable calendar time distance, this type of case ascertain- 
ment is still prone to misclassification. However, this is the PTSD 
case definition developed in September 2008 by the Department of 
Defense Interagency PTSD and TBI Standardization Committee, 
and it has been accepted by Military Health Affairs for surveillance 
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(Armed Forces Health Surveillance Center, 20 1 1 ). The DNA yields 
from the sera in this study were small, another potential limitation. 

The limitations mentioned above are offset by the rare access to 
both pre- and post-deployment biologic samples in cases and con- 
trols. While most case-control studies would not be able to infer 
whether the observed methylation patterns were a consequence of 
PTSD or whether they indicated vulnerabilities that existed among 
the cases before the onset of PTSD, our study was able to address 
both possibilities. Another strength of the study is that we uti- 
lized pyrosequencing to provide a quantitative measure of DNA 
methylation. Pyrosequencing is a quantitative real-time sequenc- 
ing technology. The light output or the Pyrosequencing signals 
(represent as Relative Light Unit or RLU) are directly proportional 
to the number of molecules being sequenced. Pyrosequencing 
QCpG software has a number of quality assessments to pass or 
to fail results for each sample being sequenced. These specifica- 
tions give information including sample quality, non-specific PCR 
amplification, non-specific sequencing, and completion of bisul- 
fite modification. Therefore, if a pyrosequencing run on a specific 
sample gave a passing result, it indicates that the DNA of that 
sample is of sufficient quality. 

Chance could account for some of the statistically significant 
findings we report, given the number of comparisons and the small 
sample size, but this was an exploratory study and the results 
should be considered as preliminary. The DoDSR, which con- 
tains over 45 million longitudinally collected serum samples of 
U.S. military service members will provide a vast resource for 
carrying out future investigations in larger populations. There is 
growing evidence that the molecular mechanisms that regulate 
DNA methylation are involved in synaptic plasticity, learning, and 
memory (Sananbenesi and Fischer, 2009). DNA is inherently sta- 
ble, compared with RNA, thus development of biomarkers based 
on DNA provide an attractive potential. Additionally, since mod- 
ifications in DNA methylation can potentially be reversed - by 
de-methylating or methylating agents - understanding the role 
of DNA methylation in PTSD has the potential to fuel novel 
approaches to PTSD treatment. 

CONCLUSION 

We evaluated DNA methylation from serum as a potential bio- 
marker for PTSD. In this exploratory study of deployed personnel, 
comparing those with PTSD and those who did not develop PTSD, 
the latter had reduced %5-mC levels of H19 and IL18 after deploy- 
ment and elevated levels of IL18 after deployment. Additionally, 
pre-deployment the people who later became cases had lower levels 



of IL18 %5-mC compared with controls. Whether these are mark- 
ers of vulnerability or resilience is a matter of speculation, but 
because of the high incidence of PTSD, it is important to study 
these patterns further. 
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